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Abstract 

For many real physico-chemical complex systems detailed mechanism includes both reversible and irreversible reac- 
tions. Such systems are typical in homogeneous combustion and heterogeneous catalytic oxidation. Most complex 
enzyme reactions include irreversible steps. The classical thermodynamics has no limit for irreversible reactions 
whereas the kinetic equations may have such a limit. We represent the systems with irreversible reactions as the limits 
of the fully reversible systems when some of the equilibrium concentrations tend to zero. The structure of the limit 
reaction system crucially depends on the relative speeds of this tendency to zero. We study the dynamics of the limit 
system and describe its limit behavior as t — > oo. The extended principle of detailed balance provides the physical 
background of this analysis. If the reversible systems obey the principle of detailed balance then the limit system with 
some irreversible reactions must satisfy two conditions: (i) the reversible part satisfies the principle of detailed balance 
and (ii) the convex hull of the stoichiometric vectors of the irreversible reactions does not intersect the linear span of 
the stoichiometric vectors of the reversible reactions. These conditions imply the existence of the global Lyapunov 
functionals and alow an algebraic description of the limit behavior. The thermodynamic theory of the irreversible 
limit of reversible reactions is illustrated by the analysis of hydrogen combustion. 

Keywords: entropy, free energy, reaction network, detailed balance, irreversibility 
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1. Introduction 



For example, let us consider a simple cycle 



1.1. The problem: non-existence of thermodynamic 
functions in the limit of irreversible reactions 

We consider a homogeneous chemical system with n 
components A,, the concentration of A, is c, > 0, the 
amount of A, in the system is N, > 0, V is the volume, 
Ni = Vc^ T is the temperature. The n dimensional vec- 
tors c = (c,) and N = (Ni) belong to the closed positive 
orthant M" in W. (IR" . (The closed positive orthant is 
the set of all vectors x e W such that x, > for all i.) 

The classical thermodynamics has no limit for irre- 
versible reactions whereas the kinetic equations have. 



k, k, h 

A\ ^ An ^ At, ^ A| 

k-i k—~) k- 3 



with the equilibrium concentrations c eq = (Cj q , c^, c^) 
and the detailed balance conditions: 

under the standard cyclic convention, here, At, + i = A\ 
and C3 + i = c\. The perfect free energy has the form 



F = V RTVci In 4r 



+ const . 
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fixed values of c^L > 0. This means that 



,, ei i 



eq 



— = -=r — > and — = — 



Let us take the fixed values of the rate constants k\ , k±2 
and A:_3. Then the limit kinetic system exists and has the 
form: 

k, fe 



k, k 2 

Ai->A 2 ^ A 3 

k-2 fc-3 



It is a routine task to write a first order kinetic equa- 
tion for this scheme. At the same time, the free energy 
function F has no limit: it tends to oo for any positive 
vector of concentrations because the term ci ln(ci/Cj q ) 
increases to oo. The free energy cannot be normalized 
by adding a constant term because the variation of the 
term c\ ln(ci/Cj q ) on an interval [0,c] with fixed c also 
increases to oo, it varies from -c^/e (for the minimizer, 
c\ = Cj q /e) to a large number c(lnc- In Cj q ) (forci = c). 

The logarithmic singularity is rather "soft" and does 
not cause a real physical problem because even for 
c^/ci = 1CT 10 the corresponding large term in the free 
energy will be just ~ 23RT per mole. Nevertheless, the 
absence of the limit causes some mathematical ques- 
tions. For example, for perfect systems with detailed 
balance under isochoric isothermal conditions the den- 
sity, 



The thermodynamic potential of a component A, can- 
not be defined in the irreversible limits when the equilib- 
rium concentration of A, tends to 0. Nevertheless, in this 
paper, we construct the universal Lyapunov functions 
for systems with some irreversible reactions. Instead 
of detailed balance we use the weaker assumption that 
these systems can be obtained from the systems with 
detailed balance when some constants tend to zero. 

1.2. The extended form of detailed balance conditions 
for systems with irreversible reactions 

Let us consider a reaction mechanism in the form of 
the system of stoichiometric equations 



a ri Ai -> YjPrjAj (r=l,...,m).. 



(2) 



/ = F/(RTV) = £ c,(ln( C/ / C ; q ) - 1) , 



(1) 



where a rl > 0, /3 r j > are the stoichiometric coeffi- 
cients. The reverse reactions with positive rate constants 
are included in the list (f2]i separately (if they exist). The 
stoichiometric vector y r of the elementary reaction is 
Jr = (jri), Jri - flri ~ a ri- We always assume that 
there exists a strictly positive conservation law, a vec- 
tor b = (hi), bi > and 2/ i>iJri - for all r. This may 
be the conservation of mass or of total number of atoms, 
for example. 

According to the generalized mass action law, the re- 
action rate for an elementary reaction (f2| is (compare to 
Eqs. (4), (7), and (14) in Ql and Eq. (4.10) in |j) 



a, = exp 



RT 



is a Lyapunov function for a system of chemical kinetics 
(here, c, is the concentration of the ;th component and 
C ; q is its equilibrium concentration for a selected value where fl . > is the activity of A,, 
of the linear conservation laws, the so-called "reference 
equilibrium"). 

This function is used for analysis of stability, exis- 
tence and uniqueness of chemical equilibrium since the 
work of Zeldovich (1938, reprinted in 1996 jH). De- 
tailed analysis of the connections between detailed bal- 
ance and the free energy function was provided in II 1911 . 
Perhaps, the first detailed proof that / is a Lyapunov 
function for chemical kinetics of perfect systems with 
detailed balance was published in 1975 12211 . 

For the irreversible systems which are obtained as 
limits of the systems with detailed balance we should 
expect the preservation of stability of the equilibrium. 
More over, one can expect existence of the Lyapunov 
functions which are as universal as the thermodynamic 
functions are. The "universality" means that these func- 
tions depend on the list of components and on the equi- 
librium concentrations but do not depend on the reaction 
rate constants directly. 



(3) 



(4) 



Here, yU, is the chemical potential and yU° is the standard 
chemical potential of the component A,. 

This law has a long history (see Jfjl 24, El 01). It 
was invented in order to meet the thermodynamic re- 
strictions on kinetics. For this purposes, according to 
the principle of detailed balance, the rate of the reverse 
reaction is defined by the same formula and its rate con- 
stant should be found from the detailed balance condi- 
tion at a given equilibrium. 

It is worth mentioning that the free energy has no 
limit when some of the reaction equilibrium constants 
tend to zero. For example, for the ideal gas the chemi- 
cal potential is ju,(c, T) = RTlnc, + /u°{T). In the irre- 
versible limit some /i9 — > oo. On the contrary, the ac- 
tivities remain finite (for the ideal gases a, = c,) and the 
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approach based on the generalized mass action law and 
the detailed balance equations w+ = w~ can be applied 
to find the irreversible limit. 

The list <(23> includes reactions with the reaction rate 
constants k r > 0. For each r we define k* = k r , w + r = w r , 
k~ is the reaction rate constant for the reverse reaction 
if it is on the list (f2]i and if it is not, w~ is the reaction 
rate for the reverse reaction if it is on the list (f2]i and if 
it is not. For a reversible reaction, K r = k^/k~ 

The principle of detailed balance for the generalized 
mass action law is: For given values k r there exists a 
positive equilibrium c£ q > with detailed balance, w* = 
w~. 

Recently, it is found the extended form of the de- 
tailed balance conditions for the systems with some ir- 
reversible reactions [12]. This extended principle of de- 
tailed balance is valid for all systems which obey the 
generalized mass action law and are the limits of the 
systems with detailed balance when some of the reac- 
tion rate constants tend to zero. It consists of two parts: 

• The algebraic condition: The principle of detailed 
balance is valid for the reversible part. (This means 
that for the set of all reversible reactions there ex- 
ists a positive equilibrium where all the elemen- 
tary reactions are equilibrated by their reverse re- 
actions.) 

• The structural condition: The convex hull of the 
stoichiometric vectors of the irreversible reactions 
has empty intersection with the linear span of the 
stoichiometric vectors of the reversible reactions. 
(Physically, this means that the irreversible reac- 
tions cannot be included in oriented cyclic path- 
ways.) 

Let us recall the formal convention: the linear span of 
empty set is {0}, the convex hull of empty set is empty. 

1.3. The structure of the paper 

In Sec. [2]we study the systems with detailed balance, 
their multiscale limits and the limit systems which sat- 
isfy the extended principle of detailed balance. The 
classical Wegscheider identities for the reaction rate 
constants are presented. Their limits when some of the 
equilibria tend to zero give the extended principle of de- 
tailed balance. 

We use the generalized mass action law for the reac- 
tion rates. For the analysis of equilibria for the general 
systems, the formulas with activities are the same as for 
the ideal systems and it is convenient to work with activ- 
ities unless we need to study dynamics. The dynamical 
variables are amounts and concentrations. In a special 



subsection l2.3l we discuss the relations between concen- 
tration and activities, formulate the main assumptions 
and present formulas for the dissipation rate. 

We introduce attractors of the systems with some ir- 
reversible reactions and study them in Sec. [3] It includes 
the central results of the paper. We fully characterize the 
faces of the positive orthant that include w-limit sets. On 
such a face, dynamics is completely degenerated (zero 
rates) or it is driven by a smaller reversible system that 
obeys classical thermodynamics. 

Hydrogen combustion is the most studied and very 
important gas reaction. It has the modest complexity: in 
the usual models there are 6-8 components and ~ 15-30 
elementary reversible reactions. Under various condi- 
tions some of these reactions are practically irreversible. 
We use this system as a benchmark in Sec.|4]and give an 
example of the correct separation of the reactions into 
reversible and irreversible part. The limit behavior of 
this system in time is described. 

In Conclusion we briefly discuss the results with fo- 
cus on the unsolved problems. 



2. Multiscale limit of a system with detailed balance 



2.1. Two classical approaches to the detailed balance 
condition 

There are two traditional approach to the description 
of the reversible systems with detailed balance. First, 
we can start from the independent rate constants of the 
elementary reactions and consider the solvability of the 
detailed balance equations as the additional condition 
on the admissible values of the rate constants. Here 
we have m constants (m should be an even number, 
m = 2t) and some equations which describe connec- 
tions between these constants. This approach was in- 
troduced by Wegscheider in 1901 12311 and developed 
further by many authors lEoi 0] ■ 

Secondly, we can select a "direct" reaction in each 
pair of mutually reverse elementary reactions. If a posi- 
tive equilibrium is known then we can find the reaction 
rate constants for the reverse reaction from the constants 
for direct reaction and the detailed balance equations. 
Therefore, the direct reaction rate constants and a set of 
the equilibrium activities form the complete description 
of the reaction. Here we have £ + n independent con- 
stants, I = m/2 rate constants of direct reactions and n 
(it is the number of components) equilibrium activities. 
For these I + n constants, the principle of detailed bal- 
ance produces no restrictions. This second approach is 
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popular in applied chemical thermodynamics and kinet- 
ics 12511 because it is convenient to work with the 
independent parameters "from scratch". 

The Wegscheider conditions appear as the necessary 
and sufficient conditions of solvability of the detailed 
balance equations. (See, for example, the textbook 
12411). Let us join the direct and reverse elementary re- 
actions and write 



when a eq , a eq 



(r=l, 



.€). 



(5) 



The stoichiometric matrix is T — (y„), y„ — - a,-i 
(gain minus loss). The stoichiometric vector y r is the 
rth row of Y with coordinates y rl = /?„■ - a r j. 

Both sides of the detailed balance equations, w* - 
w~, are positive for positive activities. The solvability of 
this system for positive activities means the solvability 
of the following system of linear equations: 

Yy ri x i = hik + r -hik; =hiK r (r =!,...€) (6) 

i 

(X{ = lna^ q ). Of course, we assume that if k* > then 
k~ > (reversibility) and the equilibrium constant K r > 
is defined for all reactions from (0. 

Proposition 1. The necessary and sufficient conditions 
for existence of the positive equilibrium a eq > with 
detailed balance is: For any solution A — (A,) of the 
system 



AT = 



i.e. ^ A r y r j = for all i 



(7) 



the Wegscheider identity holds: 



(8) 



It is sufficient to use in (0 any basis of solutions of the 
system 0: Ae {A 1 ,- - ■ ,A q }. 

2.2. Multiscale degeneration of equilibria 

Let us take a system with detailed balance and send 
some of the equilibrium activities to zero: a eq — > 
when i e I for some set of indexes /. Immediately we 
find a surprise: this assumption is not sufficient to find a 
limiting irreversible mechanism. It is necessary to take 
into account the rates of the convergency to zero of dif- 
ferent 0L q . Indeed, let us study a very simple example, 



k, k 2 
Ai ^ A 2 ^ A 3 

k-i k-2 



If a eq ,a eq — > 0, a eq /a eq = const > and a eq — 



const > then the limit system should be Aj 



A 2 
-> 0. 



A3 and we can keep ^1,-1,2 = const whereas k- 2 

If a eq ,a^ — > 0, a^/a^ — > then the limit sys- 
tem should be Ai — » A 2 — * A3 and we can keep 
ki,2 = const > whereas k-\ -2 —> 0. 

If a^ q , a eq — > 0, a eq /a eq — » then in the limit survives 
only one reaction A 2 — > A 3 (if we assume that all the 
reaction rate constants are bounded). 

We study asymptotics a eq = const x s Si , s — > 
for various values of non-negative exponents S, > 
(i = 1, . . . , «). At equilibrium, each reaction rate in the 
generalized mass action law is proportional to a power 
of e: 



w eq+ = £,tconst x e^' a ' iSi , w e r q = k r const x e 



According to the principle of detailed balance, wv q+ = 
w eq ~ and 

k + 

— = const x e (y " 6) , (9) 

K 

where 5 is the vector of exponents, 6 = (<5,-). 

There are three groups of reactions with respect to the 
given vector 5: 

l.(y r ,<5) = 0; 2.{y r ,5)<Q- 3.(y r ,S)>0. 

In the first group ((7,., 6) = 0) the ratio k+/k~ remains 
constant and we can take kf = const > 0. In the second 
group {(y r ,6) < 0) the ratio k~ /k* — > and we should 
take k~ — » whereas k^ may remain constant and pos- 
itive. In the third group ((y, , 6) > 0), the situation is 
inverse: k* /k~ — > and we can take k~ = const > 0, 
whereas k^. — > 0. 

These three groups depend on 5 but this dependence 
is piecewise constant. For every y r , three sets of 6 
are defined: (i) hyperplane (y r , 6) = 0, (ii) hemispace 
(y r , 6) < and hemispace (y r , 6) > 0. The space of 
vectors 6 is split in the subsets defined by the values of 
functions sign(y r , 6) (+1 or 0). 

We consider bounded systems, hence the negative 
values of S should be forbidden. At least one equilib- 
rium activity should not vanish. Therefore, 5j = for 
some j. Below we assume that > and 6j = for 
a non-empty set of indices Jq. Moreover, the atom bal- 
ance in equilibrium should be positive. Here, this means 
that for the set of equilibrium concentrations c eq (; e Jq) 
the corresponding values of all atomic concentrations 
are strictly positive and separated from zero. 

Let the vector of exponents, 5 = (<5,) be given and the 
three groups of reactions are found. For the reactions of 
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the third group (with (y r , 5) > 0) the direct reaction van- 
ishes in the limit s — > 0. It is convenient to transpose the 
stoichiometric equations for these reactions and swap 
the direct reactions with reverse ones. Let us perform 
this transposition. After that, a r changes over /?,., y 
transforms into —y, and the inequality (y r , 8) > trans- 
forms into (y r , 5) < 0. 

Let us summarize. We use the given vector of ex- 
ponents 6 and produce a system with some irreversible 
reactions from a system of reversible reactions and de- 
tailed balance equilibrium cL q by the following rules: 

1 . if 5, > then we assign a eq = and if S, = then 
a^ q does not change; 

2. if (y r , 6) = then kf do not change; 

3. if (y r , 6) < then we assign k~ — and k+ does 
not change; 

4. if (y r , 6) > then we assign k + r — and kr r does not 
change. (In the last case, we transpose the stoichio- 
metric equation and swap the direct reaction with 
reverse one, for convenience, y r changes to -y r and 
k~ becomes 0. Therefore, this case transforms into 
case 3.) 

This is a limit system caused by the multiscale degen- 
eration of equilibrium. The multiscale character of the 
limit a^ q = const x e 5 ' — > (for some ;) is important 
because for different values of 6 reactions may have dif- 
ferent dominant directions and the set of irreversible re- 
actions in the limit may change. 

The general form of the kinetic equations for the ho- 
mogeneous systems is 

dN 

— = V Y J w r y r , (10) 
/• 

where N, is the amount of A,, N is the vector with com- 
ponents Ni and V is the volume. 

Let us consider a limit system for the degeneration 
of equilibrium with the vector of exponents 5. For this 
system (y r , 6) < for all r and, in particular, (y r , 6) < 
for all irreversible reactions and (y r , 8) = for all 
reversible reactions. 

Proposition 2. A linear functional G#(N) = (6,N) de- 
creases along the solutions of kinetic equations for 
this limit system: dGg(N)/dt < and dGs(N)dt = 
;/ and only if all the reaction rates for the irreversible 
reactions are zero. 

Proof. Indeed, 

^fl = V^r(y r ,S)<0, (11) 

r 



because for reversible reactions (y r , 6) = 0, and for ir- 
reversible reactions w r — w+ > and (y r , S) < 0. All 
the terms in this sum are non-negative, hence it may be 
zero if and only if each summand is zero. □ 

This Lyapunov function may be used in a proof that 
the rates of all irreversible reactions in the system tend 
to with time. Indeed, if they do not tend to zero then 
on a solution of ( TTOb Gg{N(t)) —> —oo when t — > oo and 
N(t) is unbounded. Equation ( fTTT i and Proposition 
give us the possibility to prove the extended principle of 
detailed balance in the following form. Let us consider 
a reaction mechanism that includes reversible and irre- 
versible reactions. Assume that the reaction rates satisfy 
the generalized mass action law (O and the set of reac- 
tion rate constants is given. Let us ask the question: Is it 
possible to obtain this reaction mechanism and reaction 
rate constants as a limit in the multiscale degeneration 
of equilibrium from a fully reversible system with the 
classical detailed balance. The answer to this question 
gives the following theorem about the extended princi- 
ple of detailed balance. 

Theorem 1. A system can be obtained as a limit in 
the multiscale degeneration of equilibrium from a re- 
versible system with detailed balance if and only if ( i) 
the reaction rate constants of the reversible part of the 
reaction mechanism satisfy the classical principle of de- 
tailed balance and ( ii) the convex hull of the stoichio- 
metric vectors of the irreversible reactions does not in- 
tersect the linear span of the stoichiometric vectors of 
reversible reactions. 

Proof. Let the given system be a limit of a reversible 
system with detailed balance in the multiscale degener- 
ation of equilibrium with the exponent vector 5. Then 
for the reversible reactions (y r , 6) = and for the irre- 
versible reactions (y r , S) < 0. For every vector x from 
the convex hull of the stoichiometric vector of the irre- 
versible reactions (x, 6) < and for any vector y from 
the linear span of the stoichiometric vectors of the re- 
versible reactions (y, 6) = 0. Therefore, these sets do not 
intersect. The reaction rate constants for the reversible 
reactions satisfy the classical principle of detailed bal- 
ance because they do not change in the equilibrium de- 
generation and keep this property of the original fully 
reverse system with detailed balance. 

Conversely, let a system satisfy the extended princi- 
ple of detailed balance: (i) the reaction rate constants of 
the reversible part of the reaction mechanism satisfy the 
classical principle of detailed balance and (ii) the con- 
vex hull of the stoichiometric vectors of the irreversible 
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reactions does not intersect the linear span of the sto- 
ichiometric vectors of reversible reactions. Due to the 
classical theorems of the convex geometry, there exists 
a linear functional that separates this convex set from 
the linear subspace. (Strong separation of closed and 
compact convex sets.) This separating functional can be 
represented in the form (x, 9) for some vector 6. For the 
reversible reactions (y r , 6) — and for the irreversible 
reactions (y r , ff) < 0. 

It is possible to find vector 5 with this separation 
property and non-negative coordinates. Indeed, accord- 
ing to the basic assumptions, there exists a linear con- 
servation law with strongly positive coordinates. This 
is a vector b (bj > 0) with the property: (y r , b) — for 
all reactions. For any A, the vector 6 + Ab has the same 
separation property as the vector Q has. We can select 
such A that (5, = 6 t + Abj > and £, = + Ab { = for 
some i. Let us take this linear combination 5 as a vector 
of exponents. 

Let us create a fully reversible system from the ini- 
tial partially irreversible one. We do not change the re- 
versible reactions and their rate constants. Because the 
reversible reactions satisfy the classical principle of de- 
tailed balance, there exists a strongly positive vector of 
equilibrium activities a* > for the reversible reactions. 

For each irreversible reaction with the stoichiometric 
vector y r and reaction rate constant k r — k + r > we add 
a reverse reaction with the reaction rate constant 

*; = K Y](a*y y " . 

i 

For this fully reversible system the activities a* > pro- 
vide the point of detailed balance. In the multiscale de- 
generation process, the equilibrium activities depend on 
s — > as a eq = a*e s '. For the reactions with (y r , 6) — 
the reaction rate constants do not depend on s and for 
the reactions with (y,., 6) < the rate constant k~ tends 
to zero as e"™"® and k + r does not change. We return to 
the initial system of reactions in the limit s — > 0. □ 

This is a particular form of the extended principle of 
detailed balance. For more discussion see 11211 . 

2.3. Activities, concentrations and affinities 

To combine the linear Lyapunov functions Gg(N) = 
(6, N) (fTTT i with the classical thermodynamic potential 
and study the kinetic equations in the closed form we 
have to specify the relations between activities and con- 
centrations. We accept the assumption: a, = Cigi(c, T), 
where g,(c, T) > is the activity coefficient. It is a 
continuously differentiable function of c, T in the whole 



diapason of their values. In a bounded region of con- 
centrations and temperature we can always assume that 
g, > go > for some constant go. This assumption is 
valid for the non-ideal gases and for liquid solutions. It 
holds also for the "surface gas" in kinetics of hetero- 
geneous catalysis 112411 and does not hold for the solid 
reagents (see for example, analysis of carbon activity in 
the methane reforming lfl2ll ). 

The system of units should be commented. Tradi- 
tionally, a,- is assumed to be dimensionless and for per- 
fect systems = c,/c°, where c° is an arbitrary "stan- 
dard" concentration. To avoid introduction of unneces- 
sary quantities, we always assume that in the selected 
system of units, c° = 1. 

If the thermodynamic potentials exist then due to the 
thermodynamic definition of activity ©, this hypoth- 
esis is equivalent to the logarithmic singularity of the 
chemical potentials, p, = RT In c, + . . . where . . . stands 
for a continuous function of c, T (all the concentrations 
and the temperature). In this case, the free energy has 
the form 

F(N, T, V) = RT^Niilnci - 1 + f oi (c, T)) , (12) 
i 

where the functions foi(c, T) are continuously differen- 
tiable for all possible values of arguments. Functions /o, 
in the right hand side of the representation (fT2l cannot 
be restored unambiguously from the free energy func- 
tion F(N, T, V) but for a small admixture A, it is pos- 
sible to introduce the partial pressure p-, which satis- 
fies the law pi = RTcj + o(c,). This is due to the 
terms tylnc; in F. Indeed, P = -dF(N,T,V)/dV = 
RTcj + o(ci) + P L=o . Connections between the equa- 
tion of state, free energy and kinetics are discussed in 
more detail in 10, Si- 

There are several simple algebraic corollaries of the 
assumed connection between activities and concentra- 
tions. Let us consider an elementary reaction 2 a,Ai — > 
2y6,A, with ai,f3i > 0. Then, according to the general- 
ized mass action law, for any vector of concentrations c 
(fit > 0) 

1 . If, for some i, c, = then y, w(c) > 0; 

2. If, for some i, c-, = and y, < then cr, > and 
w(c) = 0. 

Similarly, for a reversible reaction 2 c/A,- ^ 

1 . If, for some i, c, = and y, > then > and 

w~(c) = 0; 

2. If, for some i, c-, = and y, < then a, > and 
w + (c) = 0. 
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These statements as well as Proposition [3] and Corol- 
lary Q] below are the consequences of the generalized 
mass action law <[3j and the connection between activ- 
ities and concentrations without any assumptions about 
extended principle of detailed balance. 

Each set of indexes J = 
the positive polyhedron, 



ij) defines a face of 



Fj — {c | Ci > for all i and c, = for i e /} . 

By definition, the relative interior of Fj, ri(Fj), consists 
of points with c, = for i e J and c,- > for i t J. 

Proposition 3. Let for a point c e ri(Fj) and an index 
i e J 

^ 7riW r {c) = . 
r 

Then this identity holds for all c e Fj. 

Proof. For convenience, let us write all the direct and 
reverse reactions separately and represent the reaction 
mechanism in the form 0. All the terms in the sum 
2,- JriWric) are non-negative, because c-, = 0. There- 
fore, if the sum is zero then all the terms are zero. The 
reaction rate w r (0 with non-zero rate constant takes 
zero value if and only if a r j > and aj = for some j. 
The equality a, = is equivalent to c, = 0. Therefore, 
w r ic) = for a point c e riiFj) if and only if there ex- 
ists j e J such that a r j > 0. If a r j > for an index j e J 
then w r ic) = for all c e Fj because cj = in Fj. □ 

We call a face Fj of the positive orthant R" invari- 
ant with respect to a set S of elementary reactions if 
IZres JrjWric) = for all c e Fj and every j e J. 

Let us consider the reaction mechanism in the form 
(0 where all the direct and reverse reactions participate 
separately. 

Corollary 1. The following statements are equivalent: 

1. Yi res JriWric) — for a point c e riiFj) and all 
indexes i 6 J; 

2. The face Fj is invariant with respect to the set of 
reactions S ; 

3. The face Fj is invariant with respect to every ele- 
mentary reaction from S ; 

4. For every r e S either y r j — for all j e J or 
a r j > Ofar some j e J. 

We aim to perform the analysis of the asymptotic be- 
havior of the kinetic equations in the multiscale degen- 
eration of equilibrium described in Sec. 12.21 For this 
purpose, we have to answer the question: how the re- 
lations between activities a,- and concentrations c,- de- 
pend on the degeneration parameter s — > 0? We do 



no try to find the maximally general appropriate answer 
to this question. For the known applications, the an- 
swer is: the relations between a,- and c, do not depend 
on e — > 0. In particular, it is trivially true for the ideal 
systems. The simple generalization, a,- = c,g,(c, T, s), 
where g,(c, T, s) > go > are continuous functions, is 
not a generalization at all, because we can use for e — > 
the limit case that does not depend on s, gjic, T) = 

gtic, r,o). 

This independence from s implies that the reversible 
part of the reaction mechanism has the thermodynamic 
Lyapunov functions like free energy. If we just delete 
the irreversible part then the classical thermodynamics 
is applicable and the thermodynamic potentials do not 
depend on e. For the generalized mass action law, the 
time derivative of the relevant thermodynamic poten- 
tials have very nice general form. Let, under given con- 
dition, the function <I>(A^, . . .) be given, where by ... is 
used for the quantities that do not change in time un- 
der these conditions. It is the thermodynamics poten- 
tial if d<$iN, . . .)/dNi = Hi. For example, it is the free 
Helmholtz energy F for V, T — const and the free Gibbs 
energy G for P, V = const. 

Let us calculate the time derivative of <f>iN, . . .) due 
to kinetic equation ( fTOb . The reaction rates are given 
by the generalized mass action law (0 with definition 
of activities through chemical potential (0. We assume 
that the principle of detailed balance holds (it should 
hold for the reversible part of the reaction mechanism 
according to the extended detailed balance conditions). 
More precisely, there exists an equilibrium with de- 
tailed balance for any temperature T, a e \T): for all r, 
w^ia eq ) = wjia eq ) = w e r q iT). It is convenient to rep- 
resent the reaction rates using these equilibrium fluxes 



+ eq 

w r = w r exp 



w r = wv q exp 



s 

i 



a ri iin-p q ) 
RT 

Priirli-rl?) 

RT 



where p q = //,(a eq , T). 

These formulas give immediately the following rep- 
resentation of the dissipation rate 



d ° _ V d ^ N ' ■ • •) AN ' 



i v dN < 



dt 



dt 



= -VRT Y (In w, + - In w;)(yv* -w~)<0. 



(13) 



The inequality holds because In is a monotone function 
and, hence, the expressions In w + r - In w~ and w + r - w~ 



7 



have always the same sign. Formulas of this kind for 
dissipation are well known since the famous Boltzmann 
//-theorem (1873 0], see also lfl3ln . The entropy in- 
crease in isolated systems has the similar form: 



dS_ 



VR y (In w+ - In w~)(w + r - w~) > . 



Let us notice that 



In w+ - lnw,. = ^ Y^^i(a ri - j3 ri ) 



RT 



The quantity -(y r ,fx) is one of the central notion of 
physical chemistry, affinity J5|. It is positive if the di- 
rect reaction prevails over reverse one and negative in 
the opposite case. It measures the energetic advantage 
of the direct reaction over the reverse one (free energy 
per mole). The activity divided by RT shows how large 
is this energetic advantage comparing to the thermal en- 
ergy. We call it the normalized affinity and use a special 
notation for this quantity: 

a,. = J*d* 

RT 



Let us apply an elementary identity 



exp a - exp b = (exp a + exp b) tanh ■ 
to the reaction rate, w r = w^-w~: 



w r = (w r + w r ) tanh — . 



a — b 



(14) 



This representation of the reaction rates gives immedi- 
ately for the dissipation rate: 



d£ 
~dt 



-VRT ^Jw + r + w~)k r tanh -y < . (15) 



In this formula, the kinetic information is collected 
in the positive factors, the sums of reaction rates 
(w* + w~), and the purely thermodynamical multipliers 
A,. tanh(A, /2) are also positive. For small |A, |, the ex- 
pression A r tanh(A, /2) behaves like A^/2 and for large 
|A r | it behaves like the absolute value, |A, |. 

So, we have two Lyapunov functions for two frag- 
ments of the reaction mechanism. For the reversible 
part, this is just a classical thermodynamic potential. 
For the irreversible part, this is a linear functional 
Gs{N) = {5,N). More precisely, the irreversible reac- 
tions decrease this functional, whereas for the reversible 
reactions it is the conservation law. Therefore, it de- 
creases monotonically in time for the whole system. 



3. Attractors 

3.1. Dynamical systems and limit points 

The kinetic equations ( TTOb do not give a complete 
representation of dynamics. The right hand side in- 
cludes the volume V and the reaction rates w r which are 
functions (0 of the concentrations c and temperature T, 
whereas in the left hand side there is N. To close this 
system, we need to express V, c and T through N and 
quantities which do not change in time. This closure 
depends on conditions. The simplest expressions ap- 
pear for isochoric isothermal conditions: V, T = const, 
c = N/V. For other classical conditions (U, V = const, 
or P, T = const, or H, P = const) we have to use the 
equations of state. There may be more sophisticated 
closures which include models or external regulators of 
the pressure and temperature, for example. 

Proposition |2 is valid for all possible closures. It is 
only important that the external flux of the chemical 
components is absent. Further on, we assume that the 
conditions are selected, the closure is done, the right 
hand side of the resulting system is continuously differ- 
entiable and there exists the positive bounded solution 
for initial data in R" and V, T remain bounded and sep- 
arated from zero. The nature of this closure is not cru- 
cial. For some important particular closures the proofs 
of existence of positive and bounded solutions are well 
known (see, for example, 12211 '). Strictly speaking, such 
a system is not a dynamical system in R" but a semi- 
dynamical one: the solutions may lose positivity and 
leave R" for negative values of time. The theory of the 
limit behavior of the semi-dynamical systems was de- 
veloped for applications to kinetic systems Jgt]. 

We aim to describe the limit behavior of the systems 
as t — > oo. Under the extended detailed balance condi- 
tion the limit behavior is rather simple and the system 
will approach steady states but to prove this behavior 
we need the more general notion of the w-limit points. 

By the definition, the w-limit points of a dynamical 
system are the limit points of the motions when time 
t — > oo. We consider a kinetic system in R" . In particu- 
lar, for each solution of the kinetic equations N(t) the set 
of the corresponding w-limit points is closed, connected 
and consists of the whole trajectories (19D, Proposition 
1.5). This means that the motion which starts from an 
w-limit point remains in R" for all time moments, both 
positive and negative. 

Proposition 4. Let N(t) be a positive solution of the ki- 
netic equation and x* be an co-limit point of this solution 
and x* — 0. then at this point i,j v > = 0. 



Proof. Let x(t) be a solution of the kinetic equations 
with the initial state x(0) - x* . All the points x{t) 
(-00 < t < 00) belong to R" . Indeed, there exists such a 
sequence tj — » 00 that iV(f/) — » jc*. For any r e (-00, 00), 
+ r) — > x(r). For sufficiently large 7, fj + t > 
and the value N(tj + r) e E". Therefore, x(r) E M" + 
(-00 < r < 00) and for any r the point x(r) is an w-limit 
point of the solution N{t). Let x* = and i,| v . = v ^ 0. 
If v > then for small |r| and r < the value of x, 
becomes negative, jc,(t) < 0. It is impossible because 
positivity. Similarly, If v < then for small r > the 
value of xi becomes negative, x,(r) < 0. It is also im- 
possible because positivity. Therefore, i,| A . =0. □ 

We use Proposition |4] in the following combination 
with Proposition[3] Let us write the reaction mechanism 
in the form (0. 

Corollary 2. If an u-limit point belongs to the relative 
interior riFj of the face Fj c R" then the face Fj is 
invariant with respect to the reaction mechanism and for 
every elementary reaction either y r j — Ofor all j e / or 
a r j > Ofor some j e /. 

Proof. If an w-limit point belongs to riFj then at this 
point all cj = for j e J due to Proposition [4] There- 
fore, we can apply Corollary Q] □ 

3.2. Steady states of irreversible reactions 

Under extended detailed balance conditions, all the 
reaction rates of the irreversible reactions are zero at 
every limit point of the kinetic equations (llOl l. due to 
Proposition [2] In this section, we give a simple combi- 
natorial description of steady states for the set of irre- 
versible reactions. This description is based on Proposi- 
tion|2]and, therefore, uses the extended detailed balance 
conditions. 

We continue to study multiscale degeneration of a 
detailed balance equilibrium. The vector of exponents 
6 = (6/) is given, > for all i and d>, = for some 
i. There are two sets of reaction. For one of them, 
{y r , 6) — and in the limit both k* > 0. In the sec- 
ond set, {y r , 6) < and in the limit we assign kf. — 
and is the same as in the initial system (before the 
equilibrium degeneration). If it is necessary, we trans- 
pose the stoichiometric equations and swap the direct 
reactions with reverse ones. 

For convenience, let us change the notations. Let 
jj be the stoichiometric vectors of reversible reactions 
with (y r , 6) = (r = 1 , . . . , h), and v/ be the stoichio- 
metric vectors for the reactions from the second set, 
(vi,6) < 0(1 = 1 ,...,s). For the reaction rates and 
constants for the first set we keep the same notations: 



w r , wf, kf.. For the second set, we use for the reaction 
rate constants qi = qt and for the reaction rates v; = v ; + . 
(They are also calculated according to the generalized 
mass action law (01.) The input and output stoichiomet- 
ric coefficients remain ay, and /?„ for the first set and for 
the second set we use the notations a], and fi v n . 

Let the rates of all the irreversible reaction be equal 
to zero. This does not mean that all the concentrations 
a, with 5\ > achieve zero. A bimolecular reaction 
A + B — » C gives us a simple example: w = ka^as 
and w = if either a a — or as — 0. On the plane 
with coordinates a a, as and with the positivity condi- 
tion, a a, as > 0, the set of zeros of w is a union of two 
semi-axes, {oa = 0, as > 0} and {oa > 0,as = 0}. 
In more general situation, the set in the activity space, 
where all the irreversible reactions have zero rates, has 
a similar structure: it is the union of some faces of the 
positive orthant. 

Let us describe the set of the steady states of the irre- 
versible reactions. Due to Proposition [2] if 2/ v/v; = 
then all v/ = 0. Let us describe the set of zeros of all vj 
in the the positive orthant of activities. 

For every / = 1 , . . . , s the set of zeros of v/ in R" is 
given by the conditions: at least for one i aZ + and 
a, = 0. It is convenient to represent this condition as a 
disjunction. Let Ji = [i \ a v n + 0}. Then the set of zeros 
of v\ an a positive orthant of activities is presented by 
the formula V/ey,( a / = 0). The set of zeros of all v; is 
represented by the following conjunction form 

A,' =I (Vfcjifa = 0)) . (16) 

To transform it into the unions of subspaces we have 
to move to a disjunction form and make some cancela- 
tions. First of all, we represent this formula as a dis- 
junction of conjunctions: 

a; =1 (v^coj = o)) 

= V, Wl ,.., <je/ , ((a„ = 0) A ... A (a it = 0)) . 

For a cortege of indexes [ii, . . . , i s ) the correspondent set 
of their values may be smaller because some values /; 
may coincide. Let this set of values be ..,,,). We can 
delete from (fTTT i a conjunction (a,, = 0) A. . . A(a, i = 0) if 
there exists a cortege {i' y , . . . , i' s ) (j!, e //) with smaller set 

of values, 5{,-, y 2 ,;). Let us check the corteges 

in some order and delete a conjunction from ( TTTb if there 
remain a term with smaller (or the same) set of index 
values in the formula. We can also substitute in (JT7J 
the corteges by their sets of values. The resulting min- 
imized formula may become shorter. Each elementary 
conjunction represents a coordinate subspace and after 
cancelations each this subspace does not belong to a 
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union of other subspaces. The final form of formula 
(Ellis 

V,(Ates,(ai = 0)), (18) 

where S j are sets of indexes, S j c { 1 , . . . , «} and for ev- 
ery two different S j, S p none of them includes another, 
Sj ^ S p . The elementary conjunction A l& $ (a; = 0) 
describes a subspace. 

The steady states of the irreversible part of the re- 
action mechanism are given by the intersection of the 
union of the coordinate subspaces ( TTSl with M". For ap- 
plications of this formula, it is important that the equal- 
ities a-, = 0, Cj = and N, = are equivalent and the 
positive orthants of the activities a,, concentrations c, or 
amounts Nj represent the same sets of physical states. 
This is also true for the faces of these orthants: Fj for 
the activities, concentrations or amounts correspond to 
the same sets of states. (The same state may corre- 
sponds to the different points of these cones, but the 
totalities of the states are the same.) 

3.3. Sets of steady states of irreversible reactions in- 
variant with respect to reversible reactions 

In this Sec. we study the possible limit behavior 
of systems which satisfy the extended detailed balance 
conditions and include some irreversible reactions. All 
the oj-limit points of such systems are steady states of 
the irreversible reactions due to Proposition|2]but not all 
these steady states may be the w-limit points of the sys- 
tem. A simple formal example gives us the couple of 
reaction: A ^ B, A + B — > C. Here, we have one re- 
versible and one irreversible reaction. The conditions of 
the extended detailed balance hold (trivially): the linear 
span of the stoichiometric vector of the reversible reac- 
tion, (-1,1,0), does not include the stoichiometric vec- 
tor of the irreversible reaction, (—1,-1,1). For the de- 
scription of the multiscale degeneration of equilibrium, 
we can take the exponents 6a = 1,6b - 1,6c -0. 

The steady states of the irreversible reaction are given 
in M" by the disjunction, (ca — 0) V (eg = 0) but only 
the points (ca = cb — 0) may be the limit points when 
t — > oo. Indeed, if ca =0 and cb > then dc^/df = 
k~CB > 0. Due to Proposition [4] this is not an w-limit 
point. Similarly, the points with ca > and eg = are 
not the w-limit points. 

Let us combine Propositions [2] [4] and Corollary [2] in 
the following statement. 

Theorem 2. Let the kinetic system satisfy the extended 
detailed balance conditions and include some irre- 
versible reactions. Then an u-limit point x* € riFj 
exists if and only if Fj consists of steady states of the 



irreversible reactions and is invariant with respect to 
all reversible reactions. 

Proof. If an w-limit point x* e riFj exists then it is a 
steady state for all irreversible reactions (due to Proposi- 
tions|2]i. Therefore, the face Fj consists of steady-states 
of the irreversible reactions (Proposition |4]i and is in- 
variant with respect to all reversible reactions (Proposi- 
tion|4]and Corollary|2]i. To prove the reverse statement, 
let us assume that Fj consists of steady states of the ir- 
reversible reactions and is invariant with respect to all 
reversible reactions. The reversible reactions which do 
not act on Cj for j e J define a semi-dynamical system 
on Fj. The positive conservation law b defines an pos- 
itively invariant polyhedron in F j. Dynamics in such a 
compact set always has w-limit points. □ 

Let us find the faces F j that contain the w-limit points 
in their relative interior riFj. According to Theorem [2] 
these faces should consist of the steady states of the ir- 
reversible reactions and should be invariant with respect 
to all reversible reactions. Let us look for the maximal 
faces with this property. For this purpose, we always 
minimize the disjunctive forms by cancelations. We do 
not list the faces that contain the w-limit points in their 
relative interior and are the proper subsets of other faces 
with this property. All the w-limit points belong to the 
union of these maximal faces. 

Let us start from the minimized disjunctive form (fT8l . 
Equation (fL8b represents the set of the steady states of 
the irreversible part of the reaction mechanism by a 
union of the coordinate subspaces A, e s,(c; = 0) in in- 
tersection with M" . It is the union of the faces, UjFs ,. If 
a face F j consists of the steady states of the irreversible 
reactions then / 2 Sj for some j. 

The following formula ( fT9] l is true on a face Fj if it 
contains w-limit points in the relative interior riFj (The- 
orem|2j: 

(c; = 0) => [(A r , rri>0 V/^x, (cj = 0)) 

A (A„ y „ <0 V jAj>0 (cj = 0))] . 

Here, c, = in Fj may be read as ; e /. Following 
the previous section, we use here the notations y n -, B r i 
and B,i for the reversible reactions and reserve v/, a\. 
and B v n for the irreversible reactions. The set of y r in 
this formula is the set of the stoichiometric vectors of 
the reversible reactions. 

The required faces F j may be constructed in an itera- 
tive procedure. First of all, let us introduce an operation 
that transforms a set of indexes S c {1, 2, . . . ,n] in a 
family of sets, &(S) = {S[, . . . , S'^. Let us take formula 



( fT9l ) and find the set where it is valid for all i e S . This 
set is described by the following formula: 

Afes |(C; = 0) A (A nrri>0 V i a >0 (Cj = 0)) 

L / J (20) 

A (A nyrt <0 V jAj>0 (Cj = 0))J . 

Let us produce a disjunctive form of this formula and 
minimize it by cancelations as it is described in Sec. 13.21 
The result is 

Vjai k (AieS^Ci = 0)) . (21) 

Because of cancelations, the sets S'j do not include one 
another. They give the result, &(S) = [S [,..., S',}. 
Each S'j e &(S ) is a superset of S , S ' 2 5. 

Let us extend the operation 6 on the sets of sets 
S = {S i, ... , S p ] with the property: 5/ <t S j for i + j. 
Let us apply & to all S , and take the union of the results: 
6()(S) = U;6(5,-). Some sets from this family may in- 
clude other sets from it. Let us organize cancelations: if 
S',S" 6 ©o(S) and 5" <z S" then retain the smallest set, 
5", and delete the largest one. We do the cancelations 
until it is possible. Let us call the final result 6(S). It 
does not depend on the order of these operations. 

Let us start from any family S and iterate the opera- 
tion 6. Then, after finite number of iterations, the se- 
quence & d (S) stabilizes: & d (S) = 6 rf+1 (S) = ... be- 
cause for any set S all sets from &(S) include S . 

The problems of propositional logic that arise in 
this and the previous section seem very similar to el- 
ementary logical puzzles yfl. In the solution we just 
use the logical distribution laws (distribution of con- 
junction over disjunction and distribution of disjunc- 
tion over conjunction), commutativity of disjunction 
and conjunction, and elementary cancelation rules like 
(A A A) o A, (AvA)e A, [A A (A V B)] o A, and 
[A V (A A B)] o A. 

Now, we are in position to describe the construction 
of all Fj that have the w-limit points on their relative 
interior and are the maximal faces with this property. 

1. Let us follow Sec. l3.2l and construct the minimized 
disjunctive form (TT~ST > for the description of the 
steady states of the irreversible reactions. 

2. Let us calculate the families of sets & d {{S j}) for 
the family of sets {S j] from ( fTST l and d = 1,2,..., 
until stabilization. 

3. Let & d ({Sj}) = e d+[ ({Sj}) = {JuJi, ■ ..J P l Then 
the family of the faces Fj. (i - 1,2, ...,p) gives 
the answer: the w-limit points are situated in riF j. 
and for each i there are w-limit points in riFj r 



3.4. Simple examples 

In this Sec, we present two simple and formal exam- 
ples of the calculations described in the previous sec- 
tions. 

l.Ai+A 2 ^A 3 +A 4 ,7 = (-1,-1,1, 1,0);A 1 +A 2 -» 
A5, v = (-1, -1,0,0, 1). The extended principle of de- 
tailed balance holds: the convex hull of the stoichiomet- 
ric vectors of the irreversible reactions consists of one 
vector j2 and it is linearly independent of j\. The in- 
put vector a for the irreversible reaction A 1 + A2 — > A 5 
is (-1,-1,0,0,0). The set J — J/ from the conjunc- 
tion form dT6b is defined by the non-zero coordinates of 
this a v : J = {1,2}. The conjunction form in this simple 
case (one irreversible reaction) loses its first conjunc- 
tion operation and is just (c\ — 0) V (C2 = 0). It is, 
at the same time, the minimized disjunction form ( fTST l 
and does not require additional transformations. This 
formula describes the steady states of the irreversible 
reaction in the positive orthant R" . For this disjunction 
form, The family of sets S = [S j) consists of two sets, 
Si = {1} andS 2 = {2}. 

Let us calculate &{S 1,2)- For both cases, ; = 1,2 there 
are no reversible reactions with y„ = 0. Therefore, one 
expression in round parentheses vanishes in (l20l . For 
S — { 1 } this formula gives 

(cj = 0) A ((c 3 = 0) V (c 4 = 0)) 

and for S = {2} it gives 

(C 2 = 0) A ((cs = 0) V (C 4 = 0)) . 

The elementary transformations give the disjunctive 
forms: 

[(ei = 0) A ((cs = 0) V (c 4 = 0))] 
o[((d = 0) A (c 3 = 0)) V (( Cl = 0) A (c 4 = 0))] , 

[(C2 = 0) A ((cs = 0) V (C 4 = 0))] 
«*[((C2 = 0) A (C 3 = 0)) V ((C 2 = 0) A (C 4 = 0))] . 

Therefore, 6(5 = {{1, 3), (1,4)}, 6(S 2 ) = 
{{2, 3}, {2,4}} and 

6({5 1 ,5 2 }) = {{1,3},{1,4},{2,3},{2,4}}. 

No cancelations are needed. The iterations of & do not 
produce new sets from {{1, 3}, {1,4}, {2, 3}, {2,4}}. In- 
deed, if ci = C3 = 0, or c\ — c\ — 0, or c 2 — C3 = 0, 
or C2 = C4 = then all the reaction rates are zero. More 
formally, for example for ©({1,3}) formula (l20l gives 

[( Cl = 0) A ((ca = 0) V (c 4 = 0))] 
A[(c 3 = 0) A (( Cl = 0) V (c 2 = 0))] . 



This formula is equivalent to (ci = 0) A (C3 = 0). There- 
fore, S({1, 3J) = {1,3}. The same result is true for {1,4}, 
{2,3}, and {2,4}. 

All the w-limit points (steady states) belong to the 
faces F {lj3 ) = {c\,c\ = c 3 = 0}, F [iA] = {c\,c\ = c 4 = 
0}, F {2 ,3) = [c\,c 2 = c 3 = 0}, or F {2 ,4} = {c\,c 2 = c 4 = 
0}. The position of the w-limit point for a solution N(t) 
depends on the initial state. More specifically, this sys- 
tem of reactions has three independent linear conserva- 
tion laws: b\ - N\ + N 2 + + Na, + 2Ns, b 2 = N\ - N 2 
and b-} = N3 -N4. For given values of these bi )2 ,3 vector 
N belongs to the 2D plane in R 5 . The intersection of 
this plane with the selected faces depends on the signs 
of b 2t y. 

• If b 2 < 0, £3 < then it intersects F{i,3] only, at one 
point N = (0, -b 2 , 0, -bj, b\+b 2 + b-}) (N5 should 
be non-negative, b\ + b 2 + b% > 0) . 

• If b 2 = 0, bj < then it intersects both F{i,3} and 
F{ 2 ,3} at one point N = (0, 0, 0, -Z?3, fri + £3) (N$ 
should be non-negative, b\ + bj > 0). 

• If b 2 < 0, £>3 = then it intersects both F(i_3} and 
F{i i 4) at one point N = (0, -b 2 , 0, 0, b\ + b 2 ) (N5 
should be non-negative, b\ + b 2 > 0). 

• If b 2 > 0, bi < then it intersects F(2,3) only, at one 
point N = (^2,0, 0, -bi,bi +b 2 + b^) (N5 should be 
non-negative, b\ + b 2 + bj, > 0). 

• If b 2 > 0, £>3 = then it intersects F^^j and F{ 2 ,4) 
at the point N = (b 2 , 0, 0, 0, b\ + b 2 ) (N$ is non- 
negative because b\ + b 2 + bi > 0). 

• lfb 2 < 0, £> 3 > then it intersects Fji^j only, at one 
point N = (0, -b 2 ,bi,0,bi + b 2 + b-}) (N 5 should be 
non-negative, b\ + b 2 + £> 3 > 0). 

• If b 2 = 0, b-i > then it intersects Fji^j and F{ 2 ,4,} 
at one point N = (0,0,& 3 ,0,&i + b^) (N 5 is non- 
negative because b\ + b^ > 0). 

• If b 2 > 0, £>3 > then it intersects F{ 2 ,4] only, at 
one point N = {b 2 , 0, bj,, 0, b\+b 2 + bj,) {N5 is non- 
negative because b\ + b 2 + bi > 0). 

As we can see, the system has exactly one w-limit point 
for any admissible combination of the values of the con- 
servation laws. These points are the listed points of in- 
tersection. 

For the second simple example, let us change the di- 
rection of the irreversible reaction. 

2. Ai + A 2 ^ A 3 + A 4 , n = (-1,-1,1,1,0), 
A 5 -> Aj + A 2 , v = (1,1,0,0,-1). The extended 



Table 1 : H 2 burning mechanism 12111 



No 


Reaction 


Stoichiometric vector 


1 


H 2 + O, ^ 20H 


(-1,-1,2,0,0,0,0,0) 


2 


H 2 + OH f± H 2 + H 


(-1,0,-1,1,1,0,0,0) 


3 


OH + O ^ O, + H 


(0,1,-1,0,1,-1,0,0) 


4 


H 2 + O f± OH + H 


(-1,0,1,0,1,-1,0,0) 


5 


O, + H +M f± H0 2 +M 


(0,-1,0,0,-1,0,1,0) 


6 


OH + H0 2 ^ 2 + H 2 


(0,1,-1,1,0,0,-1,0) 


7 


H + H0 2 f± 20H 


(0,0,2,0,-1,0,-1,0) 


8 


O + H0 2 f± 2 + OH 


(0,1,1,0,0,-1,-1,0) 


9 


20H H 2 + O 


(0,0,-2,1,0,1,0,0) 


10 


2H + M^H 2 + M 


(1,0,0,0,-2,0,0,0) 


11 


2H + H 2 ^ H 2 + H 2 


(1,0,0,0,-2,0,0,0) 


12 


2H + H 2 ^ H 2 + H 2 


(1,0,0,0,-2,0,0,0) 


13 


OH + H + M ^ H 2 + M 


(0,0,-1,1,-1,0,0,0) 


14 


H + O + M^OH + M 


(0,0,1,0,-1,-1,0,0) 


15 


20+M^0 2 + M 


(0,1,0,0,0,-2,0,0) 


16 


H + H0 2 ^ H 2 + 2 


(1,1,0,0,-1,0,-1,0) 


17 


2H0 2 ^ 2 + H 2 2 


(0,1,0,0,0,0,-2,1) 


18 


H 2 2 + M ^ 20H + M 


(0,0,2,0,0,0,0,-1) 


19 


H + H 2 2 ^ H 2 + H0 2 


(1,0,0,0,-1,0,1,-1) 


20 


OH + H 2 2 ^ H 2 + H0 2 


(0,0,-1,1,0,0,1,-1) 



principle of detailed balance holds. The steady-states 
of the irreversible reactions is given by one equation, 
C5 = 0. Formula < [2Qb gives for ©({5}) just (C5 = 0). 
The face F{5) includes w-limit points in WF(5). Dynam- 
ics on this face is defined by the fully reversible reac- 
tion system and tends to the equilibrium of the reaction 
A\ + A 2 ?± A3 + A4 under the given conservation laws. 
On this face, there exist the border equilibria, where 
c\ = c 3 = 0, or c\ = C4 = 0, or c 2 = c 3 = 0, or 
c 2 = C4 = but they are not attracting the positive solu- 
tions. 



4. Example: H2+O2 system 

For the case study, we selected the H2+O2 system. 
This is one of the main model systems of gas kinet- 
ics. The hydrogen burning gives us an example of 
the medium complexity with 8 components (Ai -Y\ 2 , 
A 2 =0 2 , A 3 =OH, A 4 =H 2 0, As =H, A 6 =0, A 7 =H0 2 , 
and A^ =H2C>2) and 2 atomic balances (H and O). For 
the example, we selected the reaction mechanism from 
1 2 111 . The literature about hydrogen burning mecha- 
nisms is huge. For recent discussion we refer to 1 3, 18 ] . 

A special symbol "M" is used for the "third body". 
It may be any molecule. The third body provides the 
energy balance. Efficiency of different molecules in this 
process is different, therefore, the "concentration" of the 
third body is a weighted sum of the concentrations of the 
components with positive weights. The third body does 
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not affect the equilibrium constants and does not change 
the zeros of the direct and inverse reaction rates but 
modifies the non-zero values of reaction rates. There- 
fore, for our analysis we can omit these terms. The 
elementary reactions 10, 11 and 12 are glued in one, 
2H^±H2, after cancelation of the third bodies, and we 
analyze the mechanism of 18 reaction. 

Under various conditions, some of the reactions are 
(almost) irreversible and some of them should be con- 
sidered as reversible. For example, let us consider the 
H2+O2 system at or near the atmospheric pressure and 
in the temperature interval 800-1200K. The reactions 
1, 2, 4, 18, 19, and 20 are supposed to be reversible 
(on the base of the reaction rate constants presented in 
|2JJ). The first question is: if these reactions are re- 
versible then which reactions may be irreversible? 

Due to the general criterion, the convex hull of the 
stoichiometric vectors of the irreversible reactions has 
empty intersection with the linear span of the stoichio- 
metric vectors of the reversible reactions. Therefore, if 
the stoichiometric vector of a reaction belongs to the lin- 
ear span of the stoichiometric vectors of the reversible 
reactions, then this reaction is reversible. Simple linear 
algebra gives that 

73,5,9 e span{yi, y 2 , 74, 7i8, 719,720} ■ 

In particular, 73 = -ji + 74, 75 = 71 - 7ig + 719, 
79 = 72 _ 74- So, the list of the reversible reactions 
should include the reactions 1, 2, 3, 4, 5, 9, 18, 19, and 
20. The reactions 6, 7, 8, 10, 11, 12, 13, 14, 15, and 17 
may be irreversible. Formally, there are 2 8 = 256 pos- 
sible combinations of the directions of these 8 reactions 
(the reactions 10, 1 1 and 12 have the same stoichiomet- 
ric vector and, in this sense, should be considered as one 
reaction). The general criterion and simple linear alge- 
bra give that there are only two admissible combinations 
of the directions of irreversible reactions: either for all 
of them k~ — or for all of them Ic} = 0. Here, the 
direct and inverse reactions and the notations kf are se- 
lected according to the TableQ] We can immediately no- 
tice that the inverse direction of all reactions is very far 
from the reality under the given conditions, for example, 
it includes the irreversible dissociation H2 — > 2H. 

Let us demonstrate in detail, how the general crite- 
rion produces this reduction from the 256 possible com- 
binations of directions of irreversible reactions to just 
2 admissible combinations. We assume that the initial 
set of reactions is spit in two: reversible reactions with 
numbers r e Jo and irreversible reactions with r e J\, 
rank(7i , 72, . . . , ye] = d, rank{7 r | r e Jo] — do. The rank 
of all vectors y r , d, must exceed the rank of the stoichio- 



metric vectors of the reversible reactions, d > do, be- 
cause if d = do then all the reactions must be reversible 
and the problem becomes trivial. 

According to 11211 . we have to perform the following 
operations with the set of stoichiometric vectors y r : 

1 . Eliminate several coordinates from all y r using lin- 
ear conservation laws. This is transfer to the inter- 
nal coordinates in span{7 r | r — !,...,{}; 

2. Eliminate coordinates from all 7,. (r e Ji) using 
the stoichiometric vectors of the reversible reac- 
tions and the Gauss-Jordan elimination procedure. 
This is the map to the quotient space span{7,. | r 

I £}/span{y r \ r e Jo}- Me denote the result as 

Jn 

3. Use the linear programming technique and analyze 
for which combinations of the signs, the convex 
hull conv{+7,. I r e J\ } does not include 0. 

In the Table [2] we present the results of the step-by- 
step elimination. First, the atomic balances give us for 
every possible stoichiometric vector 77 = (tju . . . , r]$) 
two identities: 



1. 2r]i + 773 + 2774 + Tjs + 777 + 2/78 = or 771 = 
2774 + 775 + 777 + 2t7 S ); 

2. 2772 + 773 + 774 + 776 + 2777 + 2?78 = or 772 = 
774 + 775 + 2777 + 2t7 8 ). 



-i(m 
-him 



Let us recall that the order of the coordinates (771 , . . . , 77s) 
corresponds to the following order of the components, 
(H 2 , 2 , OH, H 2 0, H, O, H0 2 , H 2 2 ). Due to these 
identities, a stoichiometric vector 77 for this mixture is 
completely defined by six coordinates (773, . . . , 77s). In 
the second column of the Table [2] these 6D vectors are 
given for all the reactions from the H2 burning mecha- 
nism (the Table [TJ. 

In five columns No. 3-7, the results of the coordinate 
eliminations are presented (and the zero-valued elimi- 
nated coordinates are omitted). Each elimination step 
may be represented as a projection: 

1 

x i-» x — Xj — 77 , 

m 

where 77, is a pivot (highlighted in bold in the column 
preceding the result of elimination), and 77 is the vec- 
tor that includes the pivot (as the zth coordinate). The 
projection operator is applied to every vector of the pre- 
vious column. At the end (the last column), all the sto- 
ichiometric vectors of the reversible reaction are trans- 
formed into zero, and the stoichiometric vectors of the 
irreversible reactions with the given direction (from the 
left to the right) are transformed into the same vector 
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Table 2: Elimination of coordinates of stoichiometric vectors for H2 burning mechanism. The reversible reactions are collected in the upper part of 
the Table. The reaction in the lower part of the table are irreversible. The group of equivalent reactions 10, 11, 12 is presented by one of them. In 
the second column, the first two coordinates (which correspond to H2 and O2) are excluded using the atomic balance. In the following columns the 
results of the coordinates elimination are presented. For each step, the pivot for elimination is underlined and highlighted in bold in the previous 
column. The eliminated coordinates at each step are named at the top of each column. Their zero values are omitted. 



No 


H 2 , 2 


OH 


H 2 o 2 


H 2 


H 


O 


1 


(2,0,0,0,0,0) 


(0,0,0,0,0) 


(0,0,0,0) 


(0,0,0) 


(0,0) 


(0) 


2 


(-1,1,1,0,0,0) 


(1,1,0,0,0) 


(1,1,0,0) 


(0,0,0) 


(0,0) 


(0) 


3 


(-1,0,1,-1,0,0) 


(0,1,-1,0,0) 


(0,1,-1,0) 


(i-1,0) 


(0,0) 


(0) 


4 


(1,0,1,-1,0,0) 


(0,1,-1,0,0) 


(0,1,-1,0) 


(1,-1,0) 


(0,0) 


(0) 


5 


(0,0,-1,0,1,0) 


(0,-1,0,1,0) 


(0,-1,0,1) 


(-1,0,1) 


(-1,1) 


(0) 


9 


(-2,1,0,1,0,0) 


(1,0,1,0,0) 


(1,0,1,0) 


(-1,1,0) 


(0,0) 


(0) 


18 


(2,0,0,0,0,-1) 


(0,0,0,0,-1) 


(0,0,0,0) 


(0,0,0) 


(0,0) 


(0) 


19 


(0,0,-1,0,1,-1) 


(0,-1,0,1,-1) 


(0,-1,0,1) 


(-1,0,1) 


(-1,1) 


(0) 


20 


(-1,1,0,0,1,-1) 


(1,0,0,1,-1) 


(1,0,0,1) 


(-1,0,1) 


(-1,1) 


(0) 


6 


(-1,1,0,0,-1,0) 


(1,0,0,-1,0) 


(1,0,0,-1) 


(-1,0,-1) 


(-1,-1) 


(-2) 


7 


(2,0,-1,0,-1,0) 


(0,-1,0,-1,0) 


(0,-1,0,-1) 


(-1,0,-1) 


(-1,-1) 


(-2) 


8 


(1,0,0,-1,-1,0) 


(0,0,-1,-1,0) 


(0,0,-1,-1) 


(0,-1,-1) 


(-1,-1) 


(-2) 


10 


(0,0,-2,0,0,0) 


(0,-2,0,0,0) 


(0,-2,0,0) 


(-2,0,0) 


(-2,0) 


(-2) 


13 


(-1,1,-1,0,0,0) 


(1,-1,0,0,0) 


(1,-1,0,0) 


(-2,0,0) 


(-2,0) 


(-2) 


14 


(1,0,-1,-1,0,0) 


(0,-1,-1,0,0) 


(0,-1,-1,0) 


(-1,-1,0) 


(-2,0) 


(-2) 


15 


(0,0,0,-2,0,0) 


(0,0,-2,0,0) 


(0,0,-2,0) 


(0,-2,0) 


(-2,0) 


(-2) 


16 


(0,0,-1,0,-1,0) 


(0,-1,0,-1,0) 


(0,-1,0,-1) 


(-1,0,-1) 


(-1,-1) 


(-2) 


17 


(0,0,0,0,-2,1) 


(0,0,0,-2,1) 


(0,0,0,-2) 


(0,0,-2) 


(0,-2) 


(-2) 



(-2). If we restore all the zeros, then the corresponding 
6D vector is (0, 0, 0, 0, -2, 0). We have to use the atomic 
balances to return to the 8D vectors. The coordinate x-j 
corresponds to HO2, x\ corresponds to H2, and xi cor- 
responds to O2, hence, 2x\ -2 = and 2x2 - 4 = 0. The 
restored 8D vector is (1 , 2, 0, 0, 0, 0, -2, 0). 

A convex combination of several copies of one vec- 
tor cannot give zero. Therefore, the structural condition 
of the extended principle of detailed balance holds. It 
holds also for the inverse direction of all the irreversible 
reactions. All other distributions of directions can pro- 
duce zero in the convex hull and are inadmissible. So, 
we have the following list of irreversible reactions that 
satisfies the extended principle of detailed balance for 
given reversible reactions. (We will not discuss the sec- 
ond list of reverse irreversible reactions because it has 
not much sense for given conditions.) 

6 OH + H0 2 -> 2 + H 2 

7 H + H0 2 -> 20H 

8 O + H0 2 -> 2 + OH 
10 2H -> H 2 

13 OH + H^H 2 

14 H + O -> OH 

15 20 -> 2 

16 H + H0 2 -> H 2 + 2 

17 2H0 2 -> 2 + H 2 2 . 



We assume that all the reaction rate constants for the 
selected directions are strictly positive. The rate of all 
these reaction vanish if and only if concentration of H, 
O and HO2 are equal to zero, cs^j = 0. Indeed, C5 = 
if and only if w\o — 0, — if and only if W15 = 0, 
a-j — if and only if wn = 0. On the other hand, all 
other reaction rates from this list are zeros if 05,6,7 = 0. 

Let us reproduce this reasoning using formulas from 
Sec. 13.21 For the Zth irreversible reaction, 7; is the set 
of indexes i for which an + 0. Let us keep for the 
irreversible reactions their numbers (6, 7, 8, 10, 13, 
14, 15, 16, 17). For them, J 6 = (3,7), J 7 = {5,7}, 
h = {6,7}, /10 = {5}, J i3 = {3,5}, J u = {5,6}, 
/i5 = {6}, /i6 = {5,7}, /i 7 = {7}. 

Formula ( TT8l gives for the steady states of the irre- 
versible reactions: 

((C3 = 0) V (C7 = 0)) A ((c 5 = 0) V (C7 = 0)) 
A((C 6 = 0) V (C7 = 0)) A (c 5 = 0) 
A((C3 = 0) V (Cj = 0)) A ((cj = 0) V (c 6 = 0)) 
A(c 6 = 0) A ((c 3 = 0) V (c 7 = 0)) A (c 7 = 0). 

It is equivalent to 

(c 5 = 0) A (c 6 = 0) A (c 7 = 0) . 

Of course, the result is the same, the face F[$fij) ( c 5,6,7 = 
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0, c, > 0) is the set of the steady states of all irreversible 
reaction. 

Let us look now on the list of reversible reactions: 

1 H 2 + 2 ^ 20H 

2 H 2 + OH ^ H 2 + H 

3 OH + O ^ 2 + H 

4 H 2 + O ^ OH + H 

5 2 + H ^ H0 2 

9 20H ^ H 2 + O 

18 H 2 2 ^ 20H 

19 H + H 2 2 ^ H 2 + H0 2 

20 OH + H 2 2 ^ H 2 + H0 2 

If the concentration OH (C3) is positive then the com- 
ponent O is produced in the reaction 9. If the concentra- 
tions of H 2 (ci) and OH (C3) both are positive then the 
component H is produced in reaction 2. If the concen- 
trations of H 2 2 (eg) and OH (C3) both are positive then 
the component H0 2 is produced in reaction 2. Due to 
the reversible reaction 18 any of two components H 2 2 
and OH produces the other component. Moreover, the 
first reaction produces OH from H 2 + 2 . This produc- 
tion stops if and only if either concentration of H 2 is 
zero (ci =0) or concentration of 2 is zero (c 2 = 0). 

This means that the set of zeros of the irreversible re- 
actions, C5 6 ; 7 = (c > 0), is not invariant with respect 
to the kinetics of the reversible reactions. This means 
that from an initial conditions on this set the kinetic tra- 
jectory will leave it unless, in addition, C3 = cs = and 
either c\ — or c 2 = 0. 

The reactions of all irreversible reactions should tend 
to zero due to Proposition [2] Therefore, the kinetic 
trajectory should approach the union of two planes, 
^1,3,5,6,7,8 = and 02,3,5,6,7,8 = (under condition c > 0). 
These planes are two-dimensional and the position of 
the state there is completely defined by the atomic bal- 
ances. 

If the concentration vector belongs to the first plane, 
then all the atoms are collected in 2 and H 2 0. It is 
possible if and only if bo > jbn- In this case, C4 = 
andc 2 = \(b - jb H ). 

If the concentration vector belongs to the second 
plane, then all the atoms are collected in H 2 and H 2 0. It 
is possible if and only if bo < \bn. In this case, c\ = bo 
and c\ = j{b H - 2b ). 

Let us reproduce this reasoning formally using the 
general formalism of Sec. 13.31 Formula [20] gives for 



©({5,6,7}) 



(C 5 


= 0)A 


k A r,y,-5>0 




>0 (Cj 


-0)) 




A 


A r , y ,. 5< 




j>o (cj 


= 0)) 


A(c 6 


= 0)A 


A,- y ,. 6> 


V;>,- 


>0 (Cj 


= 0)) 




A 


Ar,y rt «) 




,>0 (Cj 


= 0)) 


A(c 7 


= 0)A 


A r , 7r7 >o 




>0 (Cj 


= 0)) 




A 


A r>7r7 <o 


V //3 


j>o (cj 


= 0)) 



(22) 



Vectors y r in this formula participate are the sto- 
ichiometric vectors of reversible reactions (r = 
1,2, 3,4, 5, 9, 18, 19, 20). From the Table[TJwe find that 
y r5 > for r = 2, 3,4, y r5 < for r = 5, 19, y rb > for 
r = 9, y rb < for r = 3,4, y rl > for r = 5, 19,20, and 
y r i for all r. Formula (l22l transforms into 

(C 5 = 0) A ((Ci = 0)V(c 3 = 0)) A ((C3 = 0)V(C 6 = 0)) 
A(( Cl = 0)V(c 6 = 0)) A (c 7 = 0) A ((d = 0)V(c 7 = 0)) 
A(c 6 = 0) A (C3 = 0) A ((C2 = 0)V(c 5 = 0)) 
A((C 3 = 0)V(C 5 = 0)) A (C7 = 0) A ((c 2 = 0)V(c 5 = 0)) 
A((c 5 = 0)v(c 8 = 0)) A ((c 3 = 0)V(C 8 = 0)) . 

After simple transformations it becomes 

(c 3 = 0) A (c 5 = 0) A (c 6 = 0) A (c 7 = 0) . (23) 

Therefore, ©({5, 6, 7}) = {3, 5, 6, 7}. To iterate, we have 
to compute 6({3, 5, 6, 7}). For this calculation, we have 
to add one more line to formula (l22l , namely, 

A(c 3 = 0)A(A r , yr3> o V />( .. >0 (Cj = 0)) 

A (A rj7r3 <o VjfirjX) (Cj = 0)) . 

Let us take into account that y,-3 > for r = 1, 4, 18 and 
y,3 < for r = 2,3,9,20, and rewrite this formula in 
the more explicit form 

(C 3 = 0) A ((Ci = 0) V (C2 = 0)) 
A((d = 0) V (C 6 = 0)) A (eg = 0) 
A((C 4 = 0) V (C 5 = 0)) A ((C 2 = 0) V (C 5 = 0)) 
A((c 4 = 0) V (c 6 = 0)) A (c 7 = 0) . 

Let us take the conjunction of this formula with (l22l 
taken in the simplified equivalent form (l23l and trans- 
form the result to the disjunctive form. We get 

[(c 3 = 0) A (c 5 = 0) A (c 6 = 0) 
A (C 7 = 0) A (eg = 0) A (Cj = 0)] 

V[(e 3 = 0) A (e 5 = 0) A (c 6 = 0) 
A (c 7 = 0) A (c 8 = 0) A (c 2 = 0))] 



This means that 6 2 ({5,6,7}) = 6({3,5,6,7}) = 
{{1,3, 5, 6, 7, 8}, {2, 3, 5, 6, 7, 8}}. The further calcula- 
tions show that the next iteration does not change the 
result. Therefore, all the w-limit points belong to two 
faces, ^(1,3,5,6,7,81 and F{2,%5,6,7,%}- The result is the same 
as for the previous discussion. The detailed formaliza- 
tion becomes crucial for more complex systems and for 
software development. 

Let us find the vector of exponents 6 = (5,) (i = 
1, . . . , 8) from the Table [2] After all the eliminations, 
the corresponding linear functional 6 is just a value of 
the 7th coordinate: 6(x) = x-/. Its values are negative 
(-2) for all irreversible reactions and zero for all re- 
versible reactions (see the last column of the Table [2j. 
The conditions (5, y) = for the reversible reactions 
and (6, y) < for all irreversible reactions do not define 
the unique vector: if 5 satisfies these conditions then its 
linear combination with the vectors of atomic balances 
also satisfy them. Such a combination is a vector 

A5 + A H (2,0, 1,2, 1,0, l,2) + ^ o (0,2, 1, 1,0, 1,2,2) 

(25) 

under condition A > 0. This transformation of 6 does 
not change the signs of 5 on the stoichiometric vectors 
because of atomic balances. 

In our case the only coordinate remains not elimi- 
nated, x-j (the bottom part of the last column of the Ta- 
ble If, for some reaction mechanism and selected 
sets of reversible and irreversible reaction, there remain 
several {q) coordinates, then it is necessary to find q cor- 
responding functionals 6 and the space of possible vec- 
tors of exponents is (q + /(-dimensional. Here, j is the 
number of the independent linear conservation laws for 
the whole system, j = n — rankjy,}, n is the number 
of the components, \y r ) includes all the stoichiometric 
vectors for reversible and irreversible reactions. 

To find 6, we apply the elimination procedures from 
the Table|2]to an arbitrary vector y = (y,) (i = 1 , . . . , 8): 

(yi,y2,y3,y4,y5,y6,yi,y&) >-> (yi,y2,0,y4,ys,y6,yi,y&) 
"-^•(ji , 0, j 4 , J5, J6, J7, 0) h-> (yi,y 2 ,0,0,y 5 -y 4 ,y6,y 7 ,0) 
*-^(yi,y2,0,0,0,y 6 + y 5 -y 4 ,y 7 ,0) 
•— >CVi ,y~2, 0, 0, 0, 0, y 7 + y 6 + y 5 - y 4 , 0) 

(26) 

This sequence of transformations gives us the linear 
functional 

<5O0 = yi + y6 + ys - y\ ■ 

The corresponding vector of exponents 
(0,0,0,-1,1,1,1,0) should be corrected because 
its coordinates cannot be negative. Let us apply 



with A = 2 (for convenience). The coordinates 
of this combination are non-negative if and only if 
Ah > 0, Ao > and 2Au + Ao - 2 > 0. The solutions 
of these linear inequality on the (An,Ao) plane is a 
convex combination of the extreme points (corners) 
(1,0) and (0,2) plus any non-negative 2D vector: 
(iH,Ao) = ?(l,0) + (l- ? X0,2) + (#i,# 2 ),j?i,2>0and 
1 > g > 0. The corresponding vectors of exponents are 

(0, 0, 0, -2, 2, 2, 2, 0) + (g + #i)(2, 0, 1 , 2, 1 , 0, 1 , 2) 
+ (1 + # 2 )(0,4,2,2,0,2,4,4). 

At least one of the exponents should be zero. There 
are only three possibilities, S\, 62 or 64. For all other ;, 
5j > if j?i j2 > and 1 > g > 0. 

To provide any necessary atomic balance in the limit 
e — > it is necessary that two of 6, are zeros. If bo < 
4&h, then 8\ — 64 = 0. This means that #12 = 0, g = 
and 5 = (0,4,2,0,2,4,6,4). It is convenient to divide 
this 5 by 2 and write 

6 = (0,2, 2, 0,2, 2, 3, 2). 

For these exponents, the equilibrium concentrations 
tend to with the small parameter e — » (s > 0) as 



eq eq ^ eq eq 

c|7 = Cj = const, Cq = c 2 



2 eq 

e ,c QH 



eq 2 
C 3 ~ £ , 



,. et l 



eq ei 
r — r 



H 2 ~ C A ~ const ' C H _ c 5 ~ s ' c O _ "< 



„ e 'l 
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3 eq eq 
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fc ' L H : 2 L 6 



(27) 



If bo > \bn, then <S 2 = 64 — 0. This means that 
#1,2 = 0, g = 1 and 

d = (2,0,1,0,3,2,3,2). 

For these exponents, the equilibrium concentrations 
tend to with the small parameter e — > (e > 0) as 



eq eq 
c — c 
L H, L l 



2 eq 



eq eq 

c 2 = const, c QH 



eq 

c 3 ~ S, 



eq eq _ eq eq 3 eq _ eq 



eq eq 
C — C 
H0 2 7 



3 eq eq 
y r — r 
& ' H 2 2 6 



(28) 



The linear combination 2,- <WV,- decreases in time due 
to kinetic equations. This is true for any vector of expo- 
nents presented by a linear combination (EBI l (A ± 0) of 
the initial vector (0,0,0, -1, 1, 1, 1,0) with the vectors 
of the atomic balances. At the same time, any of these 
combinations give an additional linear conservation law 
for the system of reversible reactions. 

Below are several versions of this function: 
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The initial version, 6, obtained from the Table[2]is 

(6, AO = -Nh 2 o + N H + No + N H 2 ; 

Vector of exponents, calibrated by adding of the 
atomic balances (l25l l to meet the atomic balance 
conditions for bo < jba in the limit s — > is 
((J, AO = 2No 2 +2Noh+2Nh+2No+3N H o 2 +2Nh 2 o 2 ; 



• Vector of exponents, calibrated to meet the atomic 
balance conditions for bo > ^bn is, (6,N) = 
2N Hl + Nqh + 3N H + 2N + 3AW + 2AW 

All these forms differs by the combinations of the 
atomic balances d25l l and are, in this sense, equivalent. 



5. Conclusion 

The general principle of detailed balance was formu- 
lated in 1925 as follows II 1511 : "Corresponding to every 
individual process there is a reverse process, and in a 
state of equilibrium the average rate of every process is 
equal to the average rate of its reverse process." Rigor- 
ously speaking, the chemical reactions have to be con- 
sidered as reversible ones, and every step of the complex 
reaction consists of two reactions, forward and reverse 
(backward) one. However, in reality some forward or 
reverse reactions have the negligible rate. Typically, the 
complex combustion reactions, in particular, reactions 
of hydrocarbon oxidation or hydrogen combustion, in- 
clude both reversible and irreversible steps. It is a case 
in catalytic reactions as well. Although many catalytic 
reactions are globally irreversible, they always include 
some reversible steps, in particular steps of adsorption 
of gases. 

This work aims to solve the problem of the partially 
irreversible limit in chemical thermodynamics when 
some reactions become irreversible whereas some other 
reactions remain reversible. The main results in this di- 
rection are 



1 . Description of the multiscale limit of a system re- 
versible reactions when some of equilibrium con- 
centrations tend to zero ("Sec. 12.21 1. 

2. Extended principle of detailed balance for the sys- 
tems with some irreversible reactions (Theorem[T]i. 

3. The linear functional G,> that decreases in time on 
solutions of the kinetic equations under the ex- 
tended detailed balance conditions (Proposition [2] 
and Eq. (IB). 

4. The entropy production (or free energy dissipa- 
tion) formulas for the reversible part of the reac- 
tion mechanism under the extended detailed bal- 
ance conditions (Eqs. (fT3l l. (fT5ll). 



5. Description of the faces of the positive orthant 
which include the w-limit points in their relative 
interior and, therefore, description of limiting be- 
havior in time (Theorem|2]i. 

Did we solve the main problem and create the ther- 
modynamic of the systems with some irreversible re- 
action? The answer is: we solved this problem par- 
tially. We described the limit behavior but we did not 
find the global Lyapunov function that captures relax- 
ation of both reversible and irreversible parts of the sys- 
tem. The good candidate is a linear combination of the 
relevant classical thermodynamic potential and Gg but 
we did not find the coefficients. In that sense, the prob- 
lem of the limit thermodynamics remains open. 

Nevertheless, one problem is solved ultimately and 
completely: How to throw away some reverse reac- 
tions without violation of thermodynamics and micro- 
scopic reversibility? The answer is: the convex hull of 
the stoichiometric vectors of the irreversible reactions 
should not intersect with the linear span of the stoichio- 
metric vectors of the reversible reactions and the reac- 
tion rate constants of the remained reversible reactions 
should satisfy the Wegscheider identities ©. If the irre- 
versible reactions are introduced correctly then we also 
know that the closed system with this reaction mecha- 
nism goes to an equilibrium state. At this equilibrium, 
all the reaction rates are zero: the irreversible reaction 
rates vanish and the rates of the reversible reactions sat- 
isfy the principle of detailed balance. The limit equi- 
libria are situated on the faces of the positive orthant of 
concentrations and these faces are described in the pa- 
per. 

The solution of this theoretical problem is important 
for the modeling of the chemical reaction networks. 
This is because some of reactions are practically irre- 
versible. Removal of some reverse reaction from the re- 
action mechanism cannot be done independently of the 
whole structure of the reaction network; the whole reac- 
tion mechanism should be used in the decision making. 
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